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Abstract 

The Ca^"^-sensitive regulatory switch of cardiac muscle is a paradigmatic example 
of protein assemblies that communicate ligand binding through allosteric change. The 
switch is a dimeric complex of troponin C (TnC), an allosteric sensor for Ca^^, and 
troponin I (Tnl), an allosteric reporter. Time-resolved equilibrium FRET measure- 
ments suggest that the switch activates in two steps: a Tnl-independent Ca^"^-priming 
step followed by Tnl-dependent opening. To resolve the mechanistic role of Tnl in ac- 
tivation we performed stopped-flow FRET measurements of activation following rapid 
addition of a lacking component (Ca^^ or Tnl) and deactivation following rapid chela- 
tion of Ca^"^. The time-resolved measurements, stopped-flow measurements, and Ca^"^- 
titration measurements were globally analyzed in terms of a new quantitative dynamic 
model of TnC-Tnl allostery. The analysis provided a mesoscopic parameterization 
of distance changes, free energy changes, and transition rates among the accessible 
coarse-grained states of the system. The results reveal (i) the Ca^"^-induced priming 
step, which precedes opening, is the rate limiting step in activation, (ii) closing is the 
rate limiting step in deactivation, (iii) Tnl induces opening, (iv) an incompletely deac- 
tivated population when regulatory Ca^^ is not bound, which generates an accessory 
pathway of activation, and (v) incomplete activation by Ca^"^ — when regulatory Ca^"^ is 
bound, a 3:2 mixture of dynamically inter-converting open (active) and primed-closed 
(partially active) conformers is observed (15 C). Temperature-dependent stopped-flow 
FRET experiments provide a near complete thermo-kinetic parametrization of opening: 
the enthalpy change {AH = -33.4 kJ/mol), entropy change (AS' = -0.110 kJ/mol/K), 
heat capacity change {ACp = —7.6 kJ/mol/K), the enthalpy of activation (6^ = 10.6 
kJ/mol) and the effective barrier crossing attempt frequency (i^adj = 1-8 x 10^ s~^). 

Keywords allostery ; signal transduction ; energy landscape ; Markov network ; tro- 
ponin ; nonequilibrium statistical mechanics 
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Introduction 

The periodic contraction and relaxation of the heart is regulated by cytosolic Ca^^ through 
an assembly of proteins consisting, minimally, of troponin C (TnC), troponin I (Tnl), 
troponin T, tropomyosin, filamentous actin, and the head of myosin — an actin-binding 
ATPase (1-3). The TnC-Tnl assembly functions as a -sensitive regulatory switch 
as part of the larger regulatory assembly. Ca^^ binding to the single functional reg- 
ulatory site (site II) of the Ca^"^-receptor, TnC, is transduced as a delayed stochastic 
change in the isomerization state of troponin I (Tnl). The Tnl isomerization event 
relieves contractile inhibition through the release of the inhibitory region of Tnl (Tnl- 
I, residues 130-149) from actin (4). Inhibition is due the regulated association Tnl-I 
with actin that places tropomyosin in a blocking position, which inhibits inorganic 
phosphate (P^) release from myosin. release is required to form the strongly bound, 
force generating, acto- myosin complex (5). Activation of the TnC-Tnl assembly in- 
volves an inter- helical rearrangement, called opening, in the Ca^^-binding EF-hand 
(helix-loop-helix) motifs in the N-domain of TnC. The pivoting of helices B and C 
away from the central helix (D) (6) enables previously buried hydrophobic residues in 
the B, C and D helices to associate with hydrophobic residues in the regulatory helix 
of Tnl (Tnl-R, residues 150-165) that is contiguous with Tnl-I. In contrast to the fast 
skeletal muscle TnC, which has two functional regulatory Ca^"^ binding sites (sites I 
and II), in the cardiac system, Tnl-R is required for the opening of TnC (7, 8). 

The cardiac TnC-Tnl assembly is an allosteric system. Classical models of allostery, 
such as the Monod-Wyman-Changeux (9) and Koshland-Nemethy-Filmer (10) models, 
and generalized models using linked functions (11) or conditional free energies (12), 
consider systems that are in equilibrium. Allostery is the long-range coupling between 
distinct regions of a macromolecule. The classic requirements that an allosteric system 
be oligomeric and symmetric have given way to a new definition of allosteric systems — 
systems where the binding of one ligand affects the affinity of a second ligand. This 
definition includes monomeric multidomain proteins, where the second ligand is second 
protein domain, and macromolecular assemblies, where the second ligand is another 
member of the assembly. The cardiac TnC-Tnl assembly belongs to this last class. 

Over the past decade, the role of dynamics in allosteric regulation has drawn in- 
creasing interest (13). The protein is now seen as a fiuctuating entity that dynam- 
ically exchanges among a large number of microstates (13, 14). These microstates 
can be organized into regions of local stability, called conformational substates (15) or 
macrostates, that correlate with functional activity. Dynamics within a macrostate in- 
volves fiuctuating motions of individual atoms, residues, and groups of residues on the 
picosecond to microsecond timescale. These fiuctuations are the basis for the entropy 
of a macrostate (16). A second class of dynamics, which we call allosteric dynamics^ in- 
volves random jump-like transitions among macrostates on the microsecond to second 
timescale. The focus of this study is the allosteric dynamics of the cardiac TnC-Tnl 
assembly. 

Forster Resonance Energy Transfer (FRET) provides a powerful tool to study ligand 
binding-induced inter-domain distance changes (on the near A scale) in protein assem- 
blies in their native environment (17, 18). In contrast to conventional fiuorescent assays 
that report change in the local environment of the probe, FRET provides a clear metric 
for allosteric change — a change in the mean inter-probe distance. Time resolved FRET 
measurements can quantitate the distribution of inter-probe distances (19). Used in ex- 
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periments involving dynamic change through some perturbation of the system, FRET 
provides what has been termed structural kinetics (20) with structural monitoring of 
Ca^^ regulated allosteric dynamics of the TnC-Tnl assembly. 

An unfortunate property of measured perturbation-induced relaxation rates is that 
they are functions of the elementary forward and backward transition rate constants 
that govern the allosteric dynamics of a system (21). One approach to resolving the 
component elementary rate constants from observed relaxation rates is to perform a set 
of kinetic experiments that contain overlapping information. When two or more experi- 
ments jointly depend on underlying elementary transition rate parameters, a combined 
(or global) analysis of the experiments may resolve the rate parameters (22, 23). This 
resolution is usually not possible when experiments are analyzed independently and 
empirically. 

Here, we have pursued such a global strategy by using a previously character- 
ized FRET reporter system on TnC (8) in a set of experiments to parametrize the 
structural kinetics of the cardiac TnC-Tnl assembly. Global analysis of the experi- 
ments required that we formulate a mesoscopic model of TnC-Tnl allosteric dynamics 
during the activation and deactivation stages of the signaling cycle. We performed 
time resolved FRET measurements of Mg^+- or Ca^^-saturated (apo/holo) and Tnl 
bound/unbound samples; stopped-flow FRET measurements of Ca^^-induced activa- 
tion, Tnl-induced activation, and Ca^"^-chelation induced deactivation; and a FRET- 
monitored Ca^"^-titration measurement. The measurements were analyzed in terms 
of the allosteric model to provide distance changes, free energy changes, and the ele- 
mentary rate constants for transitions among the accessible coarse-grained states of the 
system. The analysis was used to differentiate between two mechanistic models of TnC- 
Tnl activation/deactivation. The experiments were used to construct the basin-limited 
portions of the free energy landscape that supports signaling. To obtain a more com- 
plete thermo-kinetic parameterization of the TnC opening transition, the stopped-flow 
FRET measurements were repeated for a set of temperatures. The recovered tran- 
sition rates for opening/closing were subjected to van't Hoff and modified Arrhenius 
analyses that quantitated the heat capacity, enthalpy change and entropy change of 
opening, the forward and backward enthalpy of activation, and the forward and back- 
ward effective barrier crossing attempt frequencies. Our measurements suggest that 
the cardiac TnC-Tnl assembly is in dynamic equilibrium among its macrostates both 
when regulatory Ca^^ is bound and unbound. When Ca2+ is bound, the probability 
that the TnC-Tnl complex is active is only 60%. 

Methods 

Sample preparation 

The construction, purification methods, and validation of TnC(F12W/N51C/C35S/C84S) 
(abbreviated name, TnC(12W/51C)) as well as tryptophan-less cardiac troponin I, 
TnI(W192F) (abbreviated name, Tnl(W-)), used in this study have been described (8). 
To remove unlabeled protein from the reaction mixture (i.e. obtain 100% labeling effi- 
ciency) the mixture was fractionated as described (24). The binary troponin complex 
was prepared by incubating TnC(12W/51C±AEDANS) (2 /iM) with an 2-fold excess 
of Tnl(W-) (4 fiM) on ice for 20 minutes. 
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Spectroscopic measurements 

Except where noted, all measurements were performed at 15 ±0.1 C. Samples were pre- 
pared in standard buffer (SB): 50 mM 3-(N-mopholino)propanesulfonic acid (MOPS), 
pH 7.0, 1 mM dithiothreitol (DTT), 5 mM MgCb, 0.2 M KCl; working buffer (WB): SB 
+ 2 mM ethylene glycol-bis-(/3-aminoethyl ether)-N,N,N',N'-tetraacetic acid (EGTA); 
or activating buffer (AB): SB + 160 fiM Ca^+ (pCa, 3.8). A stable TnC-Tnl complex 
is formed in WB with 2-fold molar excess Tnl as shown previously (25), where, us- 
ing an engineered TnI(W150), Tnl-binding-induced changes in the Trp lifetime were 
saturated with a 1.2 molar excess of TnC in the absence of Mg2+. Ca^+ titrations 
of the FRET distance were performed as described (1) with the following modifica- 
tions. A microtitrator (ISS, Champaign, IL) delivered 90 successive 5 /iL injections of 
a Ca^^-EGTA solution into reconstituted binary troponin complex ([TnC] = 1 /J.M) in 
WB (1.0 mL starting volume). The FRET donor, tryptophan, was excited at 295 nm; 
the FRET-quenched donor fluorescence was monitored at 340 nm (monochromator slit 
width, 2 fim). 

The ensemble distribution of the inter-probe distance was obtained from the time- 
resolved multi-exponential decay of the FRET donor Trpl2 as described in detail (1) 
with the following modifications. The time-resolved fluorescent decays of Trp 12 were 
collected using an IBH 5000U single photon counting time domain fluorescent lifetime 
instrument. Excitation was with a 295 nm pulsed LED. Isolated Tnl(W-) in WB 
provided background fluorescence and scattering, which was subtracted from donor- 
containing samples. (Failure to subtract background fluorescence, which is particularly 
problematic when using Trp as a FRET donor, produces reported FRET distances that 
are artificially low.) 

Perturbation-induced time-dependent FRET distances were obtained as described (26) 
using a Kintek F2004 stopped flow mixing spectrometer (1.8-ms instrument dead time). 
Deactivation kinetics: rapid mixing of a protein in AB with an equal volume of WB 
(post-mixing: 2 /xM TnC, 1 mM EGTA). Ca^^-induced activation: rapid mixing of 
the binary troponin complex in SB + minimally buffered (30 fiM EGTA was mixed 
with an equal volume of SB + 500 fiM Ca^+ (post-mixing: 2 /iM TnC, 250 fiM Ca^+). 
Tnl-induced activation: rapid mixing of a Ca^^-saturated TnC in AB with an equal 
volume of Tnl in AB. A reduced concentration of Tnl (post-mixing: 0.5 fj^M TnC, 1 
/iM Tnl) was used due to the low solubility of Tnl in 0.2 M KCl. FRET distances 
were calculated from the average of eight to ten tracings of concentration matched 
donor-only and donor- acceptor samples. Mean FRET distances recovered from the 
time-resolved data provided distance calibration standards for the stopped-flow mea- 
surements. Mock injections, where the protein was mixed with the same buffer, were 
used for calibration and to exclude the possibility of dilution-induced artifacts in the 
measured fluorescence. 

Data analysis 

The time-resolved FRET experiments were fit as described (1). Using the extensible 
global analysis software, GlobalCurve (27), the stopped flow FRET data were first in- 
dependently fit to an empirical sum of exponentials then globally fit to the kinetically 
coupled models of activation and deactivation (Eqs. 14a, 14b, and 15). Global fitting 
required custom- written routines to solve the set of coupled ordinary differential equa- 
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tions for probability distribution over the system-state space (see Modehng section) 
during activation (a) and deactivation ((i), 

^ = (aiAi+a2A2)P«(t) 

where ai is the relative species flux through path z. P = {i^\oo? ^\io? ^\ii} is the 
time-dependent marginalized probability distribution for finding TnC in state si and 
Tnl in state 52 (i.e. = (\5152)), obtained by summing over sq 

^\S1S2 ^ ^0\S1S2 + A\5lS2- (2) 

The distribution is normalized, ||^^so\5i52 II ~ "^^^ kinetic transition matrix for the 
first path of activation is 

A:_i 

Ai=| h -{k_i + k2) fc_2 |. (3) 
k2 -k-2 

The kinetic transition matrix for deactivation and the second path of activation (with 
/c^_3 = 0.1 s~^ substituted for k-s) is 



D = A2= ks -{k_3 + k2) fc-2 . (4) 






The equations were numerically integrated (Gear method, IMSL) subject to the fol- 
lowing initial conditions. The initial condition for the Ca^^-induced activation kinetics 
is the equilibrated system under deactivating conditions, P^/(0) = P|/((X)) = P|?'^ G 
A/'(D), where A/'(D) is the null space of D, obtained from the singular value decom- 
position of D, and P^?'^ is the equilibrium distribution of the protein state under 

deactivating [Ca^^]. The initial condition for Tnl-induced activation is obtained by 
setting A;2 = in Ai, giving {P^^} = {(1 + ^)"\ ^(1 + ^)"\o}. The initial 

condition for the Ca^^-chelation-induced deactivation kinetics is the equilibrium dis- 
tribution under saturating [Ca^+], P|,(0) = P^,(oo) = Pi?'"" G J\f{A). To deal with 
calibration error caused by noise in the stopped flow measurements and experimen- 
tal error in the recovered mean FRET distance from the time-resolved measurements, 
small adjustable distance correction factors were added to the models during global 
fitting. 

The observed FRET distance in the stopped flow measurements and in the Ca^"^- 
titration measurement is the population weighted average 

where = {c?\oo? ^\io? ^\ii} ^^e species-associated FRET distances. The equi- 

librium system-state distribution at the Ca^"^ chemical potential fi is given by the 
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Boltzmann equation 

where (3 = l/i?T, R is the Boltzmann gas constant, and T is absolute temperature in 
K. The set of free energies of the accessible system-states {^^^(/i)} parametrize the 
basin-limited portion of the free energy landscape. 



Modeling 

We have developed a mesoscipically-resolved theory of nonequilibrium allostery suit- 
able for describing the isomerization dynamics of TnC and Tnl as part of the cardiac 
TnC-Tnl assembly. The theory provides linked models of the activation and deactiva- 
tion components of the signaling cycle, and it can be used to combine both structural 
and kinetic information into an integrated analysis of Ca^^-mediated signaling. To 
construct a model of TnC-Tnl allostery assumptions must be made about the nature 
and number of allosteric transitions within the system. Many allosteric proteins, in- 
cluding TnC and Tnl, are adequately described as bi-meta-stable (28, 29), meaning 
they possess two macrostates: an active and an inactive isomer. The kinetics of Ca^^- 
induced structural changes in TnC-Tnl are well fit by a three-step sequential model 
involving two isomerization steps (30). The two isomerization steps can be interpreted 
as individual two-state isomerization events in the component proteins TnC and Tnl. 

Here, the allosteric protein is represented as a union of two domains — an input 
domain and an output domain. Each domain can be in one of two macrostates. Ligand 
binding to the input domain involves direct steric interaction, while the interaction 
between the input and output domains is allosteric. From the steric effects of molecular 
association (31), the state of the input domain is absolutely correlated with ligand 
binding. This assumed correlation satisfies the lock-and-key (32) and induced-fit (10) 
models. [The lock-and-key process is just an induced-fit process without structural 
change.] The system description is visualized using Scheme 1, which enumeraties all 
possible combinations of input domain, output domain, and the ligand binding status 
of the input domain. Each domain is either inactive (square) or active (circle). The 
input domain is on the left. Bound ligand is represented as a filled circle. Of the eight 
possible states, the four that are allowed are shown with a grey background. For these, 
the active input state (circle) is ligand-bound and the inactive input state (square) is 
ligand-unbound — the state of input domain is absolutely correlated with ligand binding. 
In contrast, the output domain can be either active (circle) or inactive (square) both 
when ligand is bound or unbound — it is a "pre-existing equilibrium" process (12, 33). 
Our model of allostery is thus a hybrid of pre-existing equilibrium and induced-fit 
models. The model describes allostery in all members of an assembly of allosteric 
proteins, provided that the output of a protein is viewed as the input "ligand" for the 
immediately downstream protein. In this way, inter-molecular signaling is treated as a 
steric interaction. For members of an allosteric assembly that do not function as the 
principal ligand receptor, only the output state of the protein needs to be specified 
because the state of its input is inferred from the state of the output of the upstream 
protein. 

The TnC-Tnl complex is a thermodynamic system of three components: one 
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one TnC, and one Tnl. The configuration of the entire system, called the system-state, 
is represented as a binary string, S = (so\siS2), which is a labeling scheme comprised 
of bits that specify the state of each component. The zero'th bit sq gives the ligand 
binding status of the Ca^"^-regulatory site (loop II) of TnC (0 = Ca^"^-not bound; 1 = 
Ca^"^-bound); the first bit specifies the output state of TnC, and the second (terminal) 
bit specifies the output state of Tnl. For bits representing the output state of a protein 
(bits to the right of \), denotes the inactive state and 1 denotes the active state. 

At the mesoscopic level, allosteric signaling in the TnC-Tnl assembly is a stochastic 
process because ligand-binding and protein isomerization events are random thermally- 
activated transitions. Continuous-time random processes are governed by an evolution 
equation, called the master equation, for the probability Ps{t) to find the system in 
state S at time t (34) 

= ^[W,{S'\S)Ps^{t) - W.,{S^\S)Ps{t)]. (7) 

The quantities Wp{S'\S) and W-p{S'\S) = Wp{S\S') denote, respectively the rates of 

the transitions S and S' ^ S for the elementary processes p = 1, 2, . . . , r, of which 
there are r = 8 • 7/2 = 28. Many of the fundamental properties of the nonequilibrium 
behavior of the TnC-Tnl assembly can be investigated and understood in terms of the 
Markov network of the system — a graph associated with the master equation (Eq. 7). 
In a Markov network, system-states S are represented as vertices (nodes), and the 
elementary (bi-directional) transitions p are represented as edges. 

Four carefully considered assumptions, afford a significant reduction in the number 
of elementary transitions and the number of independent kinetic/thermodynamic pa- 
rameters needed to describe Ca^^-induced signaling in the TnC-Tnl assembly, (i) Only 
one component can switch at a time. Thus, system-state changes ^ S involve tran- 
sitions in a single bit sj. Single bit transitions are represented by aj~ = {sj : ^ 1), 
and a J = {sj : 1 — > 0). The Markov network of theTnC-Tnl assembly undergoing only 
single bit transitions is shown in Fig. 2a. The number of elementary transitions are 
reduced to r = 8 • 3/2 = 12. (ii) The rate constants for each elementary transition p 
obey detailed balance 

W,{S'\S)Pl'} = W.p{S'\S)Pp, (8) 

where from Eq. 6, 

pp/pp = e-^^^^-^s'), (9) 

P^^ is the probability that the equilibrated system occupies the system-state S. We can 
identify a set of six minimal network circuits {Ci, C2, . . . Ce}, where travel along allowed 
elementary transitions p returns the system to its starting position. Each minimal 
circuit contains four nodes and consists of transitions in two bits. The cyclic flows 
and Cj^ on the circuit Ci involve, respectively, clockwise and counter-clockwise travel in 
the circuit. Conservation of free energy dictates that the net free energy change along 
each circuit Ci is zero. From conservation of free energy, Eq. 8, and Eq. 9, we find that 
for each circuit Ci the product of transition rates for must balance the product of 
transition rates for 

n = n (10) 

iec+ jec^ 
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Transition rates are thus macro scopically detail balanced. Circuits that involve transi- 
tions in all three bits can be defined, but they do not provide additional restraints on 
the transition rates. 

Constraints (i) and (ii) apply to both simple systems, such as coupled chemical 
reactions, and more complex systems, such as assemblies of allosteric proteins. The 
TnC-Tnl assembly, being a complex spatially extended system, is subject to two ad- 
ditional constraints, (iii) Each component of the assembly knows only its local envi- 
ronment. This is a defining feature of complex systems (35). It implies, for example, 
that Tnl has no direct knowledge of whether regulatory Ca^"^ is bound to TnC; rather, 
Ca^^ binding is communicated to Tnl by an allosteric change in TnC. Incomplete 
knowledge manifests as nearest-neighbor-limited influence, causing certain transition 
probabilities to be degenerate 

W{Gf\sQ . . . Sj-iSjj^i ...Sn) = W{Gf\Sj-iSjj^i), (11) 

where n = 3. Here, the rate of the transition involving bit Sj^ W{s{)\si . . . Sj-icr^Sjj^i . . . 5^), 
has been rewritten to show explicitly how the rate depends on the remaining bits Sk^j\ 
W{cr^\s{)Si . . . Sj-iSjj^i . . . Sji)- In the TnC-Tnl complex, the nearest-neighbor assump- 
tion (Eq. 11) dictates that the transition rates for Tnl switching do not directly depend 
on whether regulatory Ca^+ is bound to TnC (i.e. VI^(0\1(J^) = VF(l\l(j^)) and that 
the transition probabihties for Ca^^ binding/release do not directly depend on the 
isomerization state of Tnl (i.e. VF(cr^\siO) = W{cf^\si1)). Degenerate transition 
rates evidently originate in the constrained topology of the free energy landscape of 
the system. The forward and reverse transition probabilities Wp, W-p along each edge 
p depend on the geometry of the two-basin free energy surface of that edge (Supple- 
mental Fig. 1) as well as the basin-associated friction coefficients (36, 37). The concept 
of a one-dimensional free energy surface along the reaction coordinate of a transition is 
generalized to the composite free energy landscape^ consisting of the set of free energy 
surfaces over the reaction coordinates of all transitions p in the Markov network, (iv) 
The two system-states with inactive TnC but active Tnl — (0\01), (1\01) — are recog- 
nized as high energy states, where the Tnl-R is inserted into an unexposed hydrophobic 
pocket on TnC. Because of their high energies, the probability that these states are 
visited is negligible. These states are deleted from the network, as are the transitions 
involving them. In summary, we have introduced four assumptions that reduce the 
number of accessible network states from eight to six and constrain the transition rates 
that govern system dynamics. 

We assume that transition rates are governed by the law of mass action 

and that the transition rates obey the empirical Arrhenius Law — an exponential de- 
pendence of rate on activation energy. Out of convenience we assume that the rate of 
Ca^"^ dissociation from partially activated TnC-Tnl and fully-activated TnC-Tnl is in- 
termediate (^ 10//xs). Below, we relax this assumption and show that the assumption 
does not restrict the results. In Fig. 2b the parent Markov network (Fig. 2a) is divided 
into two sub-networks — networks for the activation and deactivation components of the 
signaling cycle. The transition rates {VFj-p} for the protein isomerizations are given in 
terms of the kinetic rate constants k±i^ i = 1, 2, 3. The activation and deactivation sub- 
networks are kinetically linked through common A:2, A:_2, ^3, and A:_3. System states 
(0\01) and (1\01) are eliminated and degenerate transition rates are shown — the tran- 
sition rates of Tnl activation ((0\10) (0\11)) and ((1\10) (1\11)) are a common 
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k2] the transition rates of Tnl deactivation ((0\10) ^ (0\11)) and ((1\10) ^ (1\11)) 
are a common k-2. W±p for the elementary Ca^^ binding transitions are not consid- 
ered. Instead, their quotients, the bi- molecular association constants Kq^ Ki^ and K2 
(in M"-"^) are defined (see Appendix B) for Ca^"^ binding/release, to the system-states 
(\00), (\10), and (\11). For both activation and deactivation, two pathways, the antic- 
ipated dominant (red) and minor (blue) pathways, are identified. Transitions common 
to both pathways are shown in green. The minor pathway for activation (deactivation) 
involves the binding (release) of Ca^"^ from a partially activated (deactivated) TnC-Tnl 
complex. The dominant and minor pathways of activation are given, respectively, by 
the kinetic schemes 

0\00 1\00 ^ 1\10 ^ 1\11 (12a) 



k-i k- 



2 



0\00 ^ 0\10 4^3, i\io ^ 1\11 (12b) 

k-3 k-2 

The dominant and minor pathways of deactivation are given, respectively, by 

0\00 ^ 0\10 ^ 0\11 4^ 1\11 (13a) 

k-3 k-2 

0\00^ 0\10 4^1^ 1\10 1\11 (13b) 

k-3 k-2 

The double arrow (<^) indicates reversible binding and release of regulatory Ca^^ by 
TnC. Activation is rightward; deactivation is leftward. 

The kinetic schemes (Eqs. 12, 13) simplify for experiments that involve rapid per- 
turbations of [Ca^^] between the extremes of activation — full saturation and full desat- 
uration. In the first path of activation (Eq. 12a), a large [Ca^^] strongly favors (1\00) 
over (0\00), and Ca^+ binding is rapid, so (0\00) can be ignored on the millisecond 
and longer time scale. In the second path of activation (Eq. 12b), rapid Cs?^ binding 
to (0\10) out-competes the deactivation of TnC through k-s- The first two steps can 
be approximated by a single transition with a forward rate constant k^ and a back- 
ward rate constant k'_^ that is small In the first path of deactivation (Eq. 13a), a 
low [Ca^^] strongly favors (0\11) over (1\11), and the step of Ca^^ release is assumed 
to not be rate limiting, so (1\11) can be ignored. In the second path of deactivation 
(Eq. 13b), a low [Ca^+] strongly favors (0\10) over (1\10). The step of Ca^+ release 
is assumed to not be rate limiting, so (1\10) can be ignored. Considering only the 
protein components {\s1S2), the schemes for the two activation pathways simplify: 

\00 ^ \10 ^ \11 (14a) 



k-i k- 



2 



\00 ^ \10 ^ \11 (14b) 

fc'_3 k-2 

The two pathways of activation will relax with different rates from differences between 
ki and ^3 and between k-i and k'_^ ~ 0. The relative flux between the two pathways 
of activation depends on the underlying rate constants for Kq and Ki as well as fci. 
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^3 and fc_3 in Eqs. 12. Relative flux can not be predicted from the reduced models 
(Eqs. 14a). 

In the two schemes for deactivation (Eq. 13), the protein system-states (\5152) and 
the kinetic transition rates are identical. The schemes are equivalent with respect 
to the measured FRET distance, which depends only on the states of the proteins. 
Considering only the protein components {\s1S2), deactivation through either (or both) 
pathway(s) will precede according to the scheme 

\00 ^ \10 ^ \11 (15) 

k-3 k-2 

The assumed ~ 10//xs transition rate for Ca^"^ release from partially-activated TnC- 
Tnl and fully-activated TnC-Tnl can be relaxed to allow faster dissociation rates. 
This leads to more complicated kinetic schemes, where Ca^^ can undergo one or more 
cycles of binding/release during activation or deactivation. In one such kinetic scheme 
of TnC-Tnl activation, Tnl activation (J2" occurs after Ca^"^ has transiently dissociated 
from TnC: 

0\00 1\00 ^ 1\10 4^ 0\10 ^ 0\11 4^ 1\11 (16) 

k-i k-2 

This scheme extends the dominant pathway of activation (Eq. 12a) by including addi- 
tional transitions of Ca^^ binding and release, but it differs from Eq. 12a because 
occurs when Ca^^ is not bound. Because the transition rates k2 and k-2 are degener- 
ate (i.e. independent of whether Ca^^ is bound to TnC) the schemes in Eq. 12a and 
Eq. 16 are kinetically equivalent with respect to the isomerization dynamics of TnC- 
Tnl (Eq. 14a). Our current measurements can not distinguish between the activation 
schemes in Eq. 12a and Eq. 16 because their isomerization dynamics are given identi- 
cally by Eq. 14a. The same arguments apply to other complicated kinetic schemes. 

Eqs. 14 and 15 are rationally reduced models for the isomerization dynamics of 
TnC and Tnl during the activation and deactivation components of the signaling cycle 
when the TnC-Tnl assembly is driven by step changes in Ca^^ between fully saturating 
and fully de-saturating conditions. Eq. 14b represents a possible accessory pathway of 
activation that may be experimentally resolved. The signaling cycle dynamics at the 
extremes of activation are determined by six rate constants: fci, A:_i, A;2, A:_2, A^a, k-^. 
Three structural configurations of TnC-Tnl are accessible: (\00), (\10), and (\11). 



Results 

Time-resolved FRET 

We employed a previously characterized FRET reporter system TnC(12W/51C*AEDANS) (8) 
to follow structural changes during the signaling cycle. The FRET assay, shown in 
Fig. 1, is sensitive to TnC opening — the inter-helical rearrangement in TnC of he- 
lices B and C relative to the central helix D. The time-resolved decays of donor-only 
TnC(12W/51C) and donor- acceptor TnC(12W/51C*AEDANS) have been reported (8) 
for isolated TnC and the TnC-Tnl complex. In this study we repeated and extended 
the measurements to include the addition of Tnl to Ca^ "^-saturated TnC. The multi- 
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exponential decays (Supplemental Fig. 2) were fit to a static Gaussian distributed 
inter-probe distance model, as described (1, 8). The recovered mean r and the stan- 
dard deviation a of the inter-probe distance distribution are shown in Fig. 3. Here, 
but not in (8), background fluorescence was subtracted from the decays. Background 
subtraction resulted in small increases in f and a reduction in a of the distance distri- 
bution (cr = l.lShw) compared to previous measurements, hw is the half width at half 
maximum that was used to quantify the breadth of the distribution in our previous 
work. 

In Fig. 3, Ca^"^-induced structural changes are apparent in both isolated TnC and 
in the binary TnC-Tnl complex. Starting with the preformed Ca^^-depleted TnC-Tnl 
complex, saturation with Ca^^ caused an apparent mean 6.6 A change in the inter- 
probe distance. Ca^"^ saturation of the previously Ca^"^-depleted TnC-only sample 
caused a small 2.2 A increase in the mean FRET distance. NMR measurements of 
the N-domain of cardiac TnC have shown a subtle 12° change in the C/D interhelical 
angle (7). The observed 2.2 A change in our FRET measurement, where the acceptor 
C51*AEDANS is in the short loop that connects the B and C hehces — the B/C linker — 
is evidently sensitive this subtle structural change. Other possibilities for the apparent 
2.2 A inter-probe distance change, such as partial opening of TnC, or altered fluo- 
rophore mobility, can not be ruled out. Subsequent addition of Tnl generated a larger 
4.4 A increase in the mean FRET distance. The sequential addition of Ca^"^ and Tnl 
produced a distance change (6.6 A) of equal magnitude to the Ca^"^-induced change 
in the binary complex, although absolute distances were 0.6 A less than for the pre- 
formed complex. These results demonstrate that activation is commutative — the order 
of Ca^"^ and Tnl addition is not important. In the analysis that follows, we take 28.9 
A as the apparent mean inter-probe distance in the fully activated binary complex. 
The standard deviation of the inter-probe distance distribution for the apo TnC-Tnl 
complex is narrower than the other samples (by ~ 1 A). 

Figure 4 shows two mechanistic models of activation that can account for the small 
change in distance observed after Ca^"^ is added to isolated TnC and the large change 
in distance that is observed after the addition of Ca^^ to the binary TnC-Tnl com- 
plex. The mechanistic models provide different structural interpretations of the al- 
losteric transitions of TnC and Tnl in the model of system dynamics (Fig. 2). In the 
population-shift-stabilization model (Fig. 4a), when Tnl is not present, Ca^"^ binding 
to TnC fails to appreciably activate TnC because the transition is energetically unfa- 
vorable, /Ci = ki/k-i < 1. Here, TnC activation is equated with TnC opening; that 
is, exposure of the hydrophobic pocket. Binding of Tnl-R to the exposed hydrophobic 
pocket stabilizes an otherwise energetically unfavorable TnC opening event. Without 
Tnl, Ca2+-bound TnC remains predominantly deactivated in the (1\00) state. Viewed 
kinetically, upon Tnl addition most Ca^^-bound TnC-Tnl must first proceed through 
the rate limiting transition ki before becoming stabilized through the k2 transition 
that out-competes the back reaction governed by ^2 > fc-i. Combining these 

results, we find that both Ca^^-induced and Tnl-induced activation are rate-limited 
by ki: k2 > k-i > ki. If TnC-Tnl switching is a population-shift-stabilization process, 
then the rate of Ca^^-induced activation and the rate of Tnl-induced activation will be 
equal The equality applies to the slowest relaxation rate of a multi-exponential decay. 

The induced-fit model interprets differently the small change in distance that is 
observed after adding Ca^^ to isolated TnC compared to the large change in distance 
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that is observed after adding Ca^ to the binary TnC-Tnl complex. In this model, 
Ca^^ binding directly and favorably activates TnC, /Ci > 1, but the TnC activation 
event, (\00) — > (\10), does not involve appreciable structural change (opening). Be- 
cause of the large /Ci, Ca^ "^-binding strongly shifts the distribution of states towards 
an energeticaly favorable "Ca^"^-primed" state of TnC (\10) that we call the primed- 
closed state. Tnl subsequently induces structural opening as Tnl-R migrates between 
the B/C and D helices of TnC, which forces TnC to open. With TnC predominantly 
in the primed-closed state, activation flux in Tnl-induced activation (Tnl addition to 
Ca^"^-bound TnC) is rate-limited only by the opening step that is governed by ^2- The 
kinetic rate of a sequential reaction cannot exceed the kinetic rate of any component 
reaction. So, the rate of Tnl-induced activation must be greater than or equal to the 
rate of Ca^"^-induced activation, which contains the isomerization steps of both TnC 
and Tnl. For the observed Ca^ "^-induced and Tnl-induced distance changes (Fig. 3), 
the induced-fit mechanism allows the rate of Tnl-induced activation to exceed the 
rate of Ca^^-induced activation. This is not the case for populaiton-shift-stabilization 
mechanism for which both Tnl-induced and Ca^ "^-induced activation are limited by ki. 

Stopped flow kinetics 

To determine which of the two mechanistic models correctly describes activation of the 
TnC-Tnl assembly, we performed a series of stopped-flow measurements of the activa- 
tion and deactivation components of the signaling cycle. In each experiment, the mean 
time-dependent FRET distance was calculated from the transient donor probe fluo- 
rescence of independently measured concentration- matched preparations of donor-only 
and donor-acceptor samples, as described (26). Trpl2 fluorescence in mock injected 
and actual samples were converted into time-dependent mean inter-probe distance us- 
ing the mean distances that were recovered from the time-resolved FRET experiments 
(Fig. 3). 

Activation kinetics were monitored in two sets of stopped-flow measurements (Fig. 5). 
For Ca^^-induced activation (Fig. 5a), distance changes were obtained after rapidly 
mixing preformed binary TnC-Tnl in a minimally Ca^^-buffered solution (30 /rM EGTA) 
with buffer containing sufficient Ca^"^ to saturate the sample (500 /iM). Empirically, 
the FRET relaxation is a bi-exponential process consisting of a very rapid transient 
(r^~^ = 1554 s~^, 48% amplitude) followed by much slower transition (t2~^ = 19 s~^, 
52% amplitude). The reduced (unweighted) x^, Xr = X^/DOF was 2.9 x 10"^. A fit 
to a three exponential function (not shown; x% 2.6 x 10~^) did not provide a statisti- 
cally significant improvement in the residuals, F = X/^(3exp)//^/?(2exp) ~ 0-87. The 68% 
confidence level is reached when F < 0.68. 

For Tnl-induced activation (Fig. 5b) distance changes were obtained from sam- 
ples of Ca^"^-presaturated TnC that were rapidly mixed with a two-fold excess of 
isolated Tnl. Rapid single exponential kinetics were observed {t~^ = 305 s"-"^, Xr = 
7.0 X 10~^). A two exponential fit (not shown; Xr = 7.0 x 10~^) did not improve 
the fitting, F = 0.99. We find that the rate of Tnl-induced activation (r~^ = 305 
g-i^ greatly exceeds the rate of Ca^"^-induced activation {r2^ = 19 s"-"^), a rate that is 
sensitive to both Ca^"^- and Tnl-activating steps. This finding is consistent with the 
induced-fit model, which, of the two mechanistic models, is the only one that allows 
the rate of cTnl-induced activation to exceed the rate of Ca^^-induced activation. If 
the transition preceded through a population-shift-stabilization mechanism then the 
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rate of Tnl-induced activation would be rate-limited to 19 s~ . The rate of overall as- 
sociation between TnC and Tnl is apparently much faster than the rate of Tnl-induced 
activation. The observed rapid single exponential transient for Tnl-induced activation 
indicates that saturation with Ca^"^ places TnC predominantly in the primed-closed 
state (1\10) that precedes opening. From the Ca^ "^-primed state, the TnC-Tnl as- 
sembly rapidly undergoes an induced- fit opening through A;2. We conclude that the 
Ca^^-priming step (with /ci < 19 s~^) is the rate limiting step in the activation of the 
TnC-Tnl assembly. The Ca^"^-priming step is energetically favorable, /Ci ^ 1 because 
there is no slow phase in Tnl-induced activation (from a population that must undergo 
a slow ki transition prior to opening through ^2). This implies that reverse reaction 
rate of Ca^^-priming is very slow: = ki/ICi ~ s~^. 

The allosteric model of activation (Fig. 2b; Eqs. 14a, 14b) provides a straight- 
forward explanation of the two-exponential process of Ca^"^-activation. The two ex- 
ponential process is interpreted as the superposition of two populations of TnC-Tnl, 
each relaxing through one of the two distinct pathways of activation. The values for 
k^_i and k^_^ that were obtained above afford preliminary estimates of ki and ^3. For 
the first pathway of activation (Eq. 14a), the dominant eigenvalue of its rate matrix 
(Eq. 3) is 

= {kik2 + k-ik-2)/{k-i + k2). (17) 

With k-i ^ 0, this simplifies to ^ fci = 19 s~^. For the second pathway of 
activation (Eq. 14b), the dominant eigenvalue of its rate matrix (Eq. 4) is 

r-i = {ksk2 + k'_^k-2)/{k'_^ + k2). (18) 

With k'_^ ~ 0, this simplifies to ks = 1554 s"-"^. The activation models provide 

a simple expression for the observed rate for Tnl-induced activation, 

= fc2 + fc_2. (19) 

Additional information needed to resolve k2 and k-2 is provided by the kinetics of 
deactivation. 

Deactivation of the TnC-Tnl assembly was monitored following rapid mixing of a 
solution containing the Ca^^ chelator EGTA (2 mM) with a Ca^"^-saturated sample of 
TnC-Tnl (Fig. 6). The relaxation transient was empirically fit as a single exponential 
decay process (r~^ = 125 s~^, = 1-7 x 10~^). A two exponential fit (not shown) did 
not statistically improve the fitting (F = 0.76). Because the scheme for deactivation 
is kinetically linked with the scheme for Tnl-induced activation (Eqs. 14, 15), the 
information in the two experiments can be combined to resolve k2 and /c_2. The 
dominant eigenvalue of the scheme of deactivation (Eq. 15) is 

= {k3k2 + k-3k-2)/{k-s + ^2). (20) 

For the system to deactivate upon the release of Ca^"^, the system must preferen- 
tially exit the primed-closed state (\10) to the left; this implies that k-s ^ k2. The 
expression for the observed rate simplifies to 



T ^ = ]Csk2 + k-2, 



(21) 
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where /C3 = ks/k-s- Eqs. 19 and 21 explain why the rate of relaxation for deactivation 
(125 s~^) is slower than the rate of Tnl-induced activation (305 s~^). In Eq. 21, k2 is 
pre- multiplied by a /C3 < 1. /C3 must be less than unity for TnC to swich off during 
deactivation. Assuming that /C3 = 0.1 we obtain estimates for k2 (200 s~^) and k-2 
(105 s"-"^). The backward reaction rate k-2 approaches the forward reaction rate A:2, 
making induced-fit TnC opening highly reversible. Under saturating Ca^^ (15 C), 
the system is an approximately 2:1 mixture of rapidly inter-converting fully activated 
(open) and partially activated (primed-closed) species. The finding that the TnC-Tnl 
does not fully activate upon Ca^"^ binding implies that the mean FRET distance of 
the Ca^^ saturated TnC-Tnl assembly (Fig. 3) is a weighted mixture of the FRET 
distances from primed-closed and open populations. 

Global analysis of equilibrium and transient FRET dis- 
tances 

Initial estimates for the rate parameters that govern the model of TnC-Tnl allostery 
(Eqs. 14 and 15) were obtained from the observed relaxation rates in the stopped-flow 
FRET measurements and assumed equilibrium constants. The structural information 
provided by FRET was ignored. To realize the full potential of the experimental infor- 
mation, we performed a global analysis of the three stopped-flow FRET measurements. 
The main task of the global fitting was to resolve the species-associated FRET distances 
and the microscopic transition rates that cause system-state populations to evolve fol- 
lowing perturbation of the system. The results of the global fitting are shown in Figs. 5 
and 6. Recovered values of the model parameters fci, A:2, A:_2, ^3, o?\oo? d\ii are 
given in Table 1. To provide stable convergence, ks d\iQ were assigned fixed 
values. The FRET distance d\iQ was fixed at 23.9 A, the mean FRET distance of the 
Ca^"^-saturated TnC-only sample from the time-resolved measurements. Parameters 
k-i and k-3 were assigned fixed values of 1.0 s~^ and 31000 s~^ because they were, 
respectively, too slow or too fast to be resolved by the data. Values were assigned based 
on the preliminary estimates of ki and ^3, the requirement that the primed-closed state 
is favored over the closed state when Ca^^ is bound, /Ci > 1, and the requirement that 
closed state is favored over the primed-closed state when Ca^^ is not bound, /C3 < 1. 

The ability of the experimental information to specify values for the model param- 
eters is quant it ated by the recovered precision of the model parameters. The precision 
of the parameters fci, k2, k-2, ^3, and d\ii were determined by projecting the 
hyper-surface (measure of deviation between model and data as a function of parame- 
ter values) along the individual parameter axes in a series of 1-D adiabatic grid searches 
(Fig. 6, inset; Table 1). The grid search results are reported as a normalized 



where Xmin the converged minimum, x^(^) values at grid search 

points. The target Xtgt — Xmini^ + ^^{^))^ is a function of the F-statistic F{5)^ 
where S is the confidence level (taken as 68%), and a ~ 1 accounts for the loss of 
one degree of freedom in the grid search. Intersection of with unity (dashed line) 
provides the upper and lower 68% (one standard deviation) error estimates of the 
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X \^) Xmin 
Xtgt - X mm 



(22) 
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parameter. In general, parameter error is non-Gaussian, and it is preferable to express 
parameter precision using upper and lower limits that non-symetrically bracket the 68% 
confidence interval. The values of parameters /ci, k2, and d\ii are determined with 
high precision (Fig. 6, inset). The parameter ks is poorly resolved since the x curve 
for ^3 is shallow; it also has more than one local minimum. This is expected since 
ks = 1618 s~^ exceeds the resolution of the stopped flow instrument (1/1.8 ms dead 
time = 555 s~^). The value of the species associated distance d\oo depends on the 
assumed value for k-i. Since the accuracy of d\QQ is conditional, its precision was not 
determined. 

The refined values of parameters k2 (181 s~^) and k-2 (117 s~^) confirmed the 
preliminary observation that Ca^^-bound TnC-Tnl is in dynamic equilibrium between 
the primed-closed and open conformations. The micro-equilibrium constant of Tnl- 
facilitated TnC opening is rigorously determined, IC2 = k2/k-2 = 1.55 (15 C). Es- 
timates for the micro-association constants that govern Ca^^-induced TnC priming 
(isomerization of A^i = ki/k-i = 18.0 and JC^ = ks/k-s = 0.052, are less reliable 
because k-i and k-s were assigned fixed values; they are, nevertheless, consistent with 
the data. The mean observed FRET distance of Ca^"^-saturated TnC-Tnl from the 
time-resolved measurements (28.9 A) is a 1.55:1 mixture of species-associated inter- 
probe distances for the primed-closed state d\iQ (mean, 23.9 A) and the open state 
d\ii (mean, 32.2 A). The mean inter-probe distance change for TnC opening is 8.3 A. 

FRET-Ca2+ titration 

From free energy conservation, the three micro-equilibrium constants that govern the 
protein isomerizations /Ci, /C2, /C3, along with any of the three equilibrium constants 
for Ca^^ binding to closed (\00), primed-closed (\10), or open (\11) system states 
(respectively, Kq^ Ki and K2) provide a complete parametrization of the free en- 
ergies {Gsi} of the relevant system-states in the macroscopic free energy landscape 
{Si^Gsi} of the TnC-Tnl assembly. To complete the thermodynamic parametriza- 
tion of the TnC-Tnl landscape, we performed an equilibrium Ca^^ titration of the 
FRET distance (Fig. 7a). For reference, the data were empirically fit to the Hill equa- 
tion, R = Rmin + (i?max " i?min)(l + lO^^^^^^-^^^^^o))-!, The rccovcrcd parameters 
(n = 1.18, pCa^o = 5.82) are consistent with previous measurements (1, 26). Through 
the Boltzmann equation (Eq. 6), the Gs^ determine how the system will equilibrate 
among all system-states Si at a given [Ca^+]. As discussed in Appendix B, the free 
energies of the ligand-unbound species depend linearly on the Ca^^ chemical potential, 
fi = —ln{10)RBTpCa: Gq\s-^s^ = G^^^^^^ + /i. The free energies of the ligand-bound 
species are independent of /x: Gi\g-^g^ = ^i\siS2' "^^^ relative free energies of the 
ligand-bound species are calculated from /c2, k-2, /ca, and k-s using 

(3AG° = -Hki/k-i), (23) 

where f3 = 1/RbT and we assign = 0. From the rate constants in Table 1, we 

obtain f3Gi\s^s^ — (0,-2.89,-3.33). The relative free energies of the ligand-unbound 
species {f3GQ\siS2) were calculated from fci, ^2, and k-2 after assigning Go\oo ~ 
where /j.' is a chemical potential offset that specifies when Go\oo ^ ^i\oo- F^^m Eq. 23 
and the rates in Table 1, we obtain, /3Go\siS2 ~ (0,2.95,2.51) + + A fit to 
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the Ca^+-FRET titration data using Eqs 2, 5, 6 (Fig. 7a) provided the value for ^/i' 
(9.56) — the remaining parameter needed to parametrize the landscape (Fig. 7b). The 
datum is sufficient to determine the remaining micro-equilibrium constants Kq^ Ki^ 
K2 that respectively determine the affinity of Ca^"^ for the system states (\00), (\10), 
(\11) (c.f. Fig. 2). Ko = exp(/3/x') = exp(9.56) = 1.4 x 10"^ M'^ From conservation of 
free energy and considering the energy changes corresponding to JCi and /C3, we obtain 
Ki = exp(9.56+2. 95+2.89) = 4.9 x 10^ M"^. We conclude that the primed-closed state 
has Ki/Kq = 300-times higher affinity for Ca^"^ than the closed state. Due to nearest- 
neighbor-limited influence, there is no relative free energy difference between the Ca^"^- 
bound and Ca^^-unbound surfaces for the opening transition, so K2 — Ki — 4.9 x 10^ 
M~^. i^app? the net (or apparent) affinity of Ca^^ for the regulatory site on TnC in 
the TnC-Tnl assembly can be calculated from the recovered landscape {5^, Gs^} using 
Eqs. 27, 29 and 6 (see Discussion section), i^app = 10^^^° = 10^'^^ = 6.0 x io^ is a 
function of the entire landscape {Si^Gs^}- The pCa^o (5.82) from the fit to the Hill 
equation closely corresponds to the recovered pK^o (5.78). 



Temperature-dependent kinetics 

To further parametrize the kinetics and thermodynamics of the opening step in activa- 
tion, we repeated the stopped flow measurements (Figs. 5 and 6) for a range of temper- 
atures T = {3.7,8, 12, 15, 16,20} C. Tnl-induced activation was single exponential at 
all temperatures with recovered r"^ = {222,254,299,305,338,362} s"^. The kinetic 
transients were globally analyzed to obtain k2{T) and k-2{T) (Fig. 8). Noise increased 
with temperature (data not shown), and the global analysis failed to converge for 20 
C. The backward rate constant k-2 showed upward deviation from a logarithmic de- 
pendence with temperature: k-2 = {49, 61, 83, 103, 117} s"^ for T = {3.7, 8, 12, 15, 16} 
C. The forward rate constant k2 showed downward deviation from a logarithmic de- 
pendence on temperature: k2 = {175, 196, 218, 191, 181} s'^ for T = {3.7, 8, 12, 15, 16} 
C. These deviations from Arrhenius Law behavior can not be attributed solely to a 
temperature-dependent energy (or enthalpy) of activation, which would cause ^2 and 
k-2 to deviate in the same direction. 

In an Arrhenius analysis of a reversible first order process, the logarithm of the 
observed relaxation rate = ^2 + k-2 is plotted against inverse temperature /3 = 
1/RbT. The slope (see Appendix A), 

?]^ = -St + Jt^AH. (24) 

df3 k2 + k-2 

depends on the enthalpy of activation of the forward reaction 6^ and the net enthalpy 
change of opening, AH. In the derivation of Eq. 24, AH and the net entropy of 
opening AS were assumed to be temperature independent. The second term in Eq. 24, 
which arises from reaction reversibility, can produce curvature in an Arrhenius plot. 
To recover the enthalpy AH and entropy AS of the opening process, we performed a 
van't Hoff analysis (Fig. 9) of the recovered temperature-dependent k2{T) and k-2{T) 
by fitting to 

/({Ai7, A5} ; 13) = ^ = e^'/^e-^^", (25) 
k-2[P) 

where /({•};/?) indicates that parameters {•} are being fit to a function / that depends 
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on f3. The rates at 16 C were excluded from the fitting. Opening is exothermic (Ai7, 
-33.4 kJ/mol) and is balanced by a loss of entropy (AS, -0.110 kJ/mol/K). Curvature 
in the data points in Fig 9 suggested that AH and AS* undergo temperature-dependent 
change over the experimental temperature range. The data were refit to Eq. 25 using 
temperature-dependent expressions for enthalpy and entropy, 

AH = AH^' + ACpiT-To) 
AS = AS^ + ACp\n{T/To) ' 

where AH^ = AH{To)^ AS^ = A5'(ro), and the heat capacity difference between the 
open and closed-primed states ACp is assumed to be temperature- independent. Fitting 
(Fig. 9, dashed line) produced maximum likelihood estimates for AH^ (-33.1 KJ/mol), 
AS"" (-0.108 KJ/mol/K), ACp (-7.6 kJ/mol/K), and Tq (282 K). The fit suggests that 
the temperature dependence of K2 = ^2/^-2 can be attributed to a large negative ACp. 
The recovered AH^ and AS^ are within 2% of the recovered temperature-independent 
AH and AS. 

To recover and the entropy- adjusted barrier crossing attempt frequency z/adj = 
z/exp {a^ /B)j is the attempt frequency, and is the entropy of activation of the 
forward reaction) the observed temperature-dependent rates = ^2 + ^-2 for Tnl- 
induced activation (c.f. Fig. 5b), were fit to (see Appendix A) 

/({.adj,5*};/?) =r-n/3) = -adje-''^* [l + e^^^'^'/^] (26) 

using the temperature-independent AH and AS obtained from the van't Hoff analysis. 
The fitting (Fig. 8, solid line) provided values for z/adj (1-8 x 10^ s"-*^) and 6'^ (10.6 
kJ/mol). Curvature in the slope is caused by reaction reversibility, a non-zero AH, 
and shifting balance of forward ^2 and reverse k-2 reaction rates. The data are well 
fit by Eq. 26, which does not include temperature-dependent changes in AH and 
AS. Using a more complicated expression with temperature-dependent AH and AS' 
is unlikely to improve the fitting results. 

The combined results of the van't Hoff and Arrhenius analyses provide the reac- 
tion enthalpy AH (-33.4 kJ/mol), the reaction entropy AS (-0.110 kJ/mol/K), the 
enthalpy of activation of the forward reaction (10.6 kJ/mol), and the entropy- 
adjusted barrier crossing attempt frequency for the forward reaction z/adj (1-8 x 10^ 
s~^). From these values we obtain the enthalpy of activation of the reverse reaction, 
^Ij^ = ^-i- — AH (44.0 kJ/mol) and the entropy- adjusted barrier crossing attempt fre- 
quency for the reverse reaction, z^adj,R = ^exp (^crl^/R^ (1.0 x lO""^^ s"-"^). z^adj,R is a 

function of the entropy of activation for the reverse reaction, = — AS, cal- 
culated using z/adj,R — ^adj exp (— AS/i?). To Confirm the fitting and calculations, 
^2 = ^adj 6xp(— /37-'-) and k-2 = z/adj,Rexp(— /35^) were calculated from the recovered 

parameters z/adj? ^Idj R' T''"' 7r- -^^S- 8, the calculated rates (dashed lines) are plot- 
ted along with the observed k2{T), k-2{T) (squares). The recovered thermo-kinetic 
parameters for the opening/closing transition are summarized in Table 2. 
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Discussion 

A previously characterized FRET reporter system TnC(12W/51C*AEDANS) (8) was 
used to probe the structural kinetics of the cardiac TnC-Tnl assembly during rapid 
activation and deactivation. In stopped-fiow kinetic measurements, FRET provides a 
meaningful measure of inter-probe distance in addition to relaxation rates. The data 
were analyzed in terms of a new non-equilibrium mesoscopic (coarse-grained) model 
of allosteric transitions in the TnC-Tnl assembly. The model captures the functional 
dynamics of the TnC-Tnl assembly in two kinetically linked allosteric models for the 
activation and deactivation stages of the signaling cycle. 

Dynamic conformational equilibrium 

Our measurements provide structural and kinetic evidence for the existence of a dy- 
namic equilibrium of macrostates in the TnC-Tnl assembly when Ca^^ is bound and 
when Ca2+ is not bound to the regulatory site (loop II) of TnC. These results pro- 
vide mechanistic insight into the phenomenon of "incomplete myofilament activation 
by Ca^+" (1, 38). In addition, the results provide, apparently, the first reported evi- 
dence of incomplete deactivation of cardiac troponin. The structural transition from 
the primed-closed state to the open state is governed by the forward rate constant ^2 
(181±{145, 222} s-^) and backward rate constant A:_2 (117±{93, 148} s"^). Parameter 
precision ± {/e, ue} is quantified by the lower and upper estimates for the 68% confi- 
dence interval. We conclude that Ca^^-bound TnC-Tnl rapidly inter-converts between 
the primed-closed and open structural conformations, spending ~65% of the time in 
the open state (T = 15 C). This finding is consistent with previously reported energetic 
measurements of Ca^^-dependent changes in the affinity of TnC for Tnl (39) and by 
high resolution structural studies of Ca^^-saturated TnC (7, 40). The affinity con- 
stants {K^^, I^a^) of native Tnl for TnC under Mg^+-saturated and Ca^+-saturated 
conditions (400mM KCl, 20 C) have been reported by Liao et. al. (39). The ratio 
/C/ = K^^/K^^ provides an estimate of the Ca^^-sensitive component of TnC-Tnl 
interaction. Equating the Ca^^-sensitive TnC-Tnl interaction with the TnC opening 
transition, we calculate, ICi = 127.0 x 10^M-V41.7 x 10^ M'^ = 3.04 (400 mM KCl, 
20 C). This is in general agreement with the isomerization constant for TnC opening 
obtained in this study, /C2 = A:2/A:_2 = 1.55 (200 mM KCl, 15 C). 

Our measurements and analysis indicate that the TnC opening causes an 8.3 A change 
in the mean inter-probe distance. This distance change is statistically indistinguishable 
from the inter-probe distance change (9.2 A) calculated from the NMR-derived Ca^^- 
bound primed-closed structure (PDB ID code 1AP4) (7) and X-ray crystallography- 
derived Ca^^-bound open structure (PDB ID code IJIE) (40) structure (JMR, manuscript 
in preperation). The inter-probe distance was obtained by modifying the high-resolution 
primed-closed and open structures in-silico to incorporate the donor and acceptor 
FRET probes (Fig. 1) and performing all- atom molecular dynamics simulations to 
sample the conformations of the FRET probes. 

When regulatory Ca^"^ is not bound to TnC, the TnC-Tnl assembly may exist in 
the partially activated primed-closed state in addition to the closed state. Activation 
through the accessory pathway in our allosteric model (Fig. 4b (blue); Eq 12b) occurs 
through Ca^^-binding to a pre-existing population of TnC-Tnl in the primed-closed 
state. Our kinetic measurements detect activation through this pathway. Solzin et 
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al., using an environmentally sensitive fluorophore to monitor activation kinetics, also 
detected a rapid transient during Ca^^-activation (41). They interpreted this transient 
as an additional step that precedes the rate-limiting (Ca^^-priming) step in a serial 
reaction scheme for Ca^"^-activation. From the FRET-provided distance information in 
our measurements, we found that the rapid Ca^"^-induced transient (1554 s"-"^) involves 
substantial (48% of 6.6 A) distance change. This large distance change is associated 
with opening (not priming). We have shown that opening temporally follows the slow 
priming transition. The accessory pathway of activation thus bypasses the slow priming 
step in activation through Ca^"^ binding to a pre-existing population of TnC-Tnl in the 
primed-closed state. The accessory pathway of activation suggests that the TnC-Tnl 
assembly does not completely deactivate to the closed conformation when regulatory 
Ca^"^ is not bound. 

Mesoscopic modeling 

Our mesoscopic model of TnC-Tnl allostery (Fig. 2b) recognizes transition rate degen- 
eracy due to nearest-neighbor-limited influence, and it recognizes the kinetic linkages 
between activation and deactivation stages of the signaling cycle. Transition rate de- 
generacy arises from the complex nature of the system. Simple systems such as chem- 
ical networks — networks of chemical reactions of diffusing molecules — do not exhibit 
transition rate degeneracy. The model shows, for example, that the elementary tran- 
sition rate parameters that govern the accessory pathway of activation are the same 
rate parameters that govern the kinetics of deactivation. These linkages impact both 
structural and kinetic aspects of the measured inter-probe distance relaxations in the 
two activation measurements, the deactivation measurement and the Ca^"^-titration 
measurement. Recognizing these linkages enabled recovery of the model parameters 
from the measurements through a global analysis of the data because the observed 
experimental relaxation rates jointly depend upon the linked model parameters. Inde- 
pendent empirical fitting of stopped flow experiments ignores the intensity information 
in the measurements and the inherent kinetic linkages in the system. In the compart- 
mental analysis of kinetic transitients — the most common type of analysis — intensity 
information is not utilized because (i) analytical expressions for intensity are compli- 
cated functions of the transition rates and initial conditions, (ii) initial conditions must 
be known or assumed, and (iii) the species-associated intensities, including interme- 
diate species, must be known. The lack of constraint when data are independently 
fit and empirically fit to sums of exponentials makes drawing biological conclusions 
problematic. In contrast, global analysis of a diverse set of experiments in terms of a 
single model facilitates strong conclusions because global analysis is self-constraining — 
recovered model parameters must account for the observed behavior in all experiments. 

Induced-fit opening 

Two mechanistic models of TnC-Tnl allostery (Fig. 4) — the population-shift-stabilization 
model and induced-fit model — are consistent with the small Ca^"^-induced change in 
isolated TnC and the large distance change observed upon saturation of TnC-Tnl with 
Ca^^ (Fig- 3). These models provide different interpretations with respect to struc- 
ture, kinetics, and energetics of the transitions in the model of TnC-Tnl allostery 
(Fig. 2b). The population-shift-stabilization model interprets the isomerization of 
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TnC as a structural opening event that, in the absence of Tnl, is energetically un- 
favorable and statistically improbable. Tnl switching, through Tnl-R association with 
the opening-induced hydrophobic pocket in TnC, stabilizes the otherwise unfavorable 
transition. The induced-fit model interprets the TnC isomerization as a structurally 
subtle, but energetically favorable, transition that generates a Ca^"^-primed state with 
a structure that is closed. TnC opening is the Tnl-mediated opening of permissive 
Ca^^-primed TnC. The two models have different kinetic signatures. Of the two, only 
the induced-fit model allows the rate of Tnl-induced activation to exceed the rate 
of Ca^"^-induced activation. The mechanism is diagnosed by an observed rate of Tnl- 
induced activation (305 s~^, Fig 5b) that greatly exceeds the slow rate of Ca^"^-induced 
activation (19 s~^. Fig 5a). 

Theories of allosterism assume that ligand agonists directly stabilize the active 
state of their receptor (reviewed in (42)). In the population-shift-stabilization model, 
Ca^^ binding does not directly stabilize the receptor TnC. Instead TnC is indirectly 
stabilized by Ca^^ in a process that relies on a free energy drop provided by Tnl. In 
the induced-fit mechanism, Ca^"^ binding directly stabilizes the primed-closed confor- 
mation of TnC, which we find has a 300-times greater affinity for Ca^"^ than the closed 
conformation. Our finding that the cardiac TnC-Tnl assembly switches as an induced- 
fit process questions the belief that energetics must correlate with structural change, 
i.e. small energy drops produce little structural change; large energy drops produce 
appreciable movement. TnC opening, which involves an 8.3 A inter-probe distance 
change, is, thermodynamically, marginally favorable (AG = —1.1 kJ/mol, 15 C). In 
contrast, the Ca^^-priming step, which produces a subtle 2.2 A inter-probe distance 
change, is thermodynamically very favorable (AG = —6.9 kJ/mol, 15 C). Mechanisti- 
cally, induced-fit opening appears to involve the replacement of intra-TnC hydrophobic 
interactions with inter-protein TnC-Tnl hydrophobic interactions as the Tnl-R diffu- 
sively creeps deeper into the TnC hydrophobic pocket. TnC opening is coupled to, and 
limited by, Tnl isomerization. 

Free energy landscape 

The mesoscopic free energy landscape consists of the set of one dimensional free energy 
surfaces along the reaction coordinates of all system-state transitions (edges in the 
Markov network, Fig. 2a). System-states (nodes in the Markov network) are regions 
of local stability (low energy wells) in the free energy landscape. Collectively, the 
system-states Si and the macroscopic free energies of the system states Gs^ define the 
macroscopic free energy landscape^ or free energy landscape {Si^Gsi} of the system. 
The free energy landscape of the TnC-Tnl assembly (excluding the high energy states 
(\01)) is shown in Fig 7b. The free energy landscape is of fundamental interest because 
it provides the weighting scheme needed to calculate the macroscopic measure of any 
observable property of the equilibrated system. Through the Boltzmann equation 
(Eq. 6), the Gs^ determine the equilibrium probability P^^(/i) that the system occupies 
system-state Si under an externally imposed Ca^^ chemical potential /i. The ensemble 
measurement of any observable O^. of the equilibrated system is the inner-product of 
the observable with the equilibrium probability distribution 




(27) 



{S^} 
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The O-based pCa^o = lii(lO) is obtained from the where 

\ Si / max \ 'S'z / min 

The FRET distance d^. (defined in Methods) in a Ca^^ titration experiment (c./. 
Fig. 7a) is one such system observable O^.; many others exist. The distance-^^Caso or 
pd5o(5.78) was obtained by solving Eq. 28 for the observable (i^.. It is generally assumed 
that the sensitivity to Ca^^ of some physical observable of myofilament regulation (e.g. 
force, probe fluorescence) is equal to the net affinity of the TnC regulatory site for 
Ca^"^. To test this assumption we examined the correspondence between the pd^o and 
the net Ca^"^-binding-p(7a5o or pK^Q. The Ca^ "^-binding status B5. of system-state Si 
is defined 

(S5'.)(/i) was calculated using Eqs. 27, 28, and 29. Significantly, the traces of ((i5'.)(/i) 

and (^^^(/i) superimpose, and the recovered ^^50 and pK^o are identically 5.78 (data 
not shown). This supports the assumption that the sensitivity to Ca^"^ of a system 
observable is a measure of net Ca^^ affinity for the regulatory site of TnC. 

In addition to its quantitative value, the free energy landscape (Fig 7b) provides a 
visual representation of how activation/deactivation occurs. The Boltzmann equation 
(Eq. 6) dictates that the system will respond to any perterbation by preferentially 
populating the regions of lowest free energy In the TnC-Tnl assembly, as [Ca^^] is 
raised, the G^. of the Ca^^-unbound system states (0\5i52) (green plane) rise above 
the Ca^^-bound system states (l\5i52) (blue plane), which are stationary. Population 
density (represented as sphere radius in Fig 7b) is greatest for those states with lowest 
free energy at a given [Ca^^] (pCa). Activity is defined as a change in the probability 
that Tnl, the reporting component of the assembly, is either active (§2 = 1) or inactive 
(52 = 0). Visualized in terms of the free energy landscape, activation is rightward 
movement of density towards the Tnl-active states (0\5il) and (l\5il). Deactivation 
is leftward movement of density away from the Tnl-active states. At low [Ca^^], the 
system is mostly deactivated in (0\00) with very little partially activated species (0\10). 
At high [Ca^^] a 3:2 mixture of primed-closed (l\5iO) and open (l\5il) system-states 
occur. 

Although the micro-equilibrium constants /Ci = ki/k-i = 18.0 and /C3 = ks/k-^ = 
0.052 were determined inexactly, we can confirm that their values jointly satisfy the re- 
quirement that the work done on the free energy landscape by Ca^^ binding not exceed 
the input of free energy to the assembly from Ca^"^ binding. From the Ca^"^ titration 
measurement , Ca^+ binding provides AGin = ln{10)RTpCa^o = 31.9 kJ/mol of free 
energy into the TnC-Tnl assembly (15 C). This free energy is the maximum amount of 
energy available to change the free energy landscape of the assembly. Kinetic linkages 
governing the Tnl isomerization, (\10) ^ (\11), dictate that the bound and unbound 
surfaces have the same AG for the Tnl isomerization, i.e. Gi^n — Gi\io = Go\ii — Go\io- 
However, the Ca^^-bound and Ca^^-unbound surfaces of the free energy landscape have 
different AG for the TnC isomerization, (\00) ^ (\10). Indeed, without some differ- 
ence in AG between the Ca^"^-bound and Ca^"^-unbound surfaces, no Ca^"^-induced 
population shift would occur for the protein states {\s1S2). From the estimated /Ci 
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and /C3, the total relative change between the unbound and bound free energy surfaces 
is AGout = —RTln{ICs/ICi) = 14.0 kJ/mol. The efficiency of coupling the free energy 
of ligand binding into a change in the free energy landscape, e = AGout/ AGin — 0.44. 
As required, e < 1. 

Thermo-kinetics of TnC-Tnl opening 

From temperature-dependent stopped flow kinetics, we have determined the enthalpy 
AH (-33.4 kJ/mol), entropy AS (-0.110 kJ/mol/K), and heat capacity change ACp 
(-7.6 kJ/mol/K) of the isolated opening transition in the activation of the cardiac 
TnC-Tnl assembly (Table 2). The measured values are for the opening step in ac- 
tivation, which is distinct from other processes such as Ca^"^-binding, Ca^"^-priming, 
and overall TnC-Tnl association. Opening is promoted by a favorable decrease in en- 
thalpy that is partially offset by a loss of entropy. Like cardiac TnC, calmodulin is 
a member of the super-family of bi-lobed Ca^^ binding proteins with two EF-hand 
domains per lobe. The entropy change for the binding of a peptide of myosin light- 
chain kinase to calmodulin has been studied by deuterium NMR relaxation (43) , where 
side chain order parameters provide an upper limit on AS (-0.4752 kJ/mol/K). Our 
measured AS' (-0.110 kJ/mol/K) for Tnl-R-induced TnC opening follows this trend 
of entropy loss upon complexation. Our measured AH and AS' are opposite those 
observed for tension generating step in psoas muscle, which is endothermic and dis- 
ordering (44, 45). Large negative ACp have been consistently observed for ligand 
binding- induced isomerizations of proteins. Examples include ACp = —3.5 kJ/mol/K 
for DNA binding by the transcription factor Oct-1 (46) and ACp = —3.1 kJ/mol/K for 
NAD+ binding to glyceraldehyde-3-phosphate dehydrogenase (47). Large negative Cp 
indicate sequence-specific recognition (48). In most cases, ACp are measured as average 
quantities for coupled equilibria involving combinations of bound/unbound ligand and 
active/inactive protein. In this context, ACp is properly regarded as an apparent heat 
capacity changeACp^^ with contributions from the intrinsic ACp between the inactive 
and active protein states and the ligand-induced redistribution of the macrostate pop- 
ulation (49). Our measurements have resolved the true ACp for the opening transition 
of the TnC-Tnl assembly. 

The measured Ai7, AS', and ACp can be used to recover the distribution of enthalpy 
(energy) in the open and primed-closed macrostates. The enthalpy distribution arises 
from structural fluctuations — inter-converting microstates — within each macrostate. 
Each macrostate consists of a set of N microstates mi, 7712, . . . , mA^, and each microstate 
has an enthalpy H{mi). We define a set of discrete enthalpy intervals i^i, 7^2, • • • , ^^oo 
of constant width A. We cluster the microstate enthalpies H{mi) into the enthalpy 
levels Hi] each level contains some number gi microstates. The probability of finding 
the macrostate in the enthalpy interval Hi is p{Hi) — Qi/N . The continuous probability 
distribution of enthalpy is obtained in the limit, p{H)dH = Hira^o p{Hi). Evidently, 
p{H) is the probability that a macrostate has enthalpy in the interval [H, H + dH) 
given that H is available (^ e~^^, a decreasing function in H) times the probability 
that H is available (ri(i7), an increasing function in H). p{H) is therefore a peaked 
function. The enthalpy distribution p[H) is parameterized by its moments 
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Ignoring moments in p{H) greater than two, p{H) is Gaussian 

{H - {H)Y 



P{H) « ^= exp 



2ct2 



(30) 



with mean {H) — iJ^bs variance related to the enthalpy fluctuations 5H by cr^ = 
{SH'^) = {H'^) — {H)'^ . The variance is related to the heat capacity through 

= RT'^Cp. (31) 

By introducing the subscript {k) into p{H)^ a, etc., we emphasize that the 

quantities are properties of the macrostate {k) — they are mesoscopic quantities — not 
macroscopic quantities of the system as a whole. The enthalpy distribution of 
macrostate {k) is thus parametrized by = {H)(j^^ and Cp^j^y 

We now show how the measured AS' for opening is used to resolve the heat ca- 
pacities, hence the width of the enthalpy distributions, of the open and primed-closed 
macrostates. The macrostate, a mesoscopic entity, is a sufficiently large collection of 
microstates that thermodynamic quantities can be appropriately defined for it by sta- 
tistical thermodynamics (50). As with macroscopic entropy, the mesoscopic entropy 
S'(/.) of macrostate (k) is a functional of the probability distribution of enthalpy, 

/oo 
-OO 

For Gaussian distributed p(/.)(i7) (Eq. 30) with standard deviation (7(/.), we obtain (51), 

S^k)=R\nV2^(T^k)' (32) 

Quite remarkably, (j(^) is completely determined by S'(^). From Eqs. 31 and 32, the 
entropy change for opening, AS = 5'(5q\ii) — ^(soVio)? ^ simple expression of the heat 
capacities of the open and primed-closed macrostates. 



AS = kslnJ^^. (33) 
The specific heat capacities of the open and primed-closed states are 

^P{so\10) = ^2AS/R_ I ' 

Evaluating, we obtain Cp(^^^\^^Q^^ = 7.6 kJ/mol and C'p^^^^-,^-,^^J ~ lO"-"^^ kJ/mol. The 
extremely low value of C^pi^gQ^n) is due to the highly negative AS, causing e^AS/i? _^ q 
The enthalpy distribution of the primed-closed state P(^so\io) {H) is broad with cr(5Q\io) = 
72 kJ/mol (T = 15 C). Compared to p(^g^\iQ^{H), the enthalpy distribution of the open 
state P(so\ii){H) is almost delta-correlated ct^^^^iq) ^ kJ/mol, and shifted to lower 
enthalpy (by -33.4 kJ/mol). 

We have obtained, additionally, the enthalpy of activation for opening 5^ (10.6 
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kJ/mol) and the entropy-adjusted, or effective, barrier crossing attempt frequency for 
opening z/adj = ^exp (cr^/i?) (1.8 x 10^ s"-*^), where is the entropy of activation, z/adj is 
very slow, and the rate of TnC opening is hmited by a large loss of entropy needed 
to reach the transition state. We conclude that TnC-Tnl opening involves a large 
transient loss (a^) of configurational freedom and a net loss (AS') of configurational 
freedom. 

Finally, we note that the third law of thermodynamics demands that Cp be non- 
negative (50). From Eq. 34b we obtain the constraint 

> 0, (35) 



^2AS/R _ I 

which is satisfied by our results. This provides a useful check for experimental results. 



Conclusions 

Our results indicate the presence of dynamic exchange on two levels: dynamic exchange 
among macrostates of the cardiac TnC-Tnl regulatory complex (the closed, primed- 
closed, and open macrostates) and dynamic exchange among microstates within each 
macrostate. A statistical description is necessary to explain the macroscopic (phys- 
iologic) behavior of cardiac regulation because activation and deactivation are never 
complete. The chance that the system occupies an active or inactive macrostate is 
only more or less probable. One thus speaks of the cardiac regulatory cycle in terms 
of shifting macrostate populations determined by a macroscopic free energy landscape 
that is modulated by the cytosolic concentration of Ca^+. 
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Here we derive expressions that allow an experimental thermo-kinetic parametrization 
of the energy landscape of a bi-meta-stable system that reversibly and stochastically 
jumps between two states, A and B (supplemental Fig. 1), according to the kinetic 
scheme 

A ^ B. (36) 

k-i 

The reversible transition is of interest not only in single step reactions but also as a 
component step of multi-step reactions in coupled allosteric transitions. An Arrhenius 
analysis of temperature-dependent kinetic data, where the logarithm of the observed 
reaction rate is plotted against inverse temperature, is typically limited to the analysis 
of an irreversible reaction 

A^BorA^^'B (37) 

either because, ki ^ A;_i, or, ki <^ A:_i, or because separate ki and k-i can be 
resolved through external information (45). We wish to perform an Arrhenius analysis 
of the reversible first order process in Eq. 36. A general expression is derived for the 
dependence of the observed relaxation rate, = ki + on inverse temperature. 
The expression applies to both reversible and irreversible reactions. 

We assume that both forward and reverse isomerizations in Eq. 36 are governed by 
the van't Hoff- Arrhenius law 

kj^i = Uj^i exp(-/3AG]_.) (38) 

with effective barrier crossing attempt rate uj^i from Kramers' reaction rate theory 
(36, 53, 54), barrier crossing free energy AG^^^, and /? = 1/RT. Barrier crossing free 
energies are defined by 

7^ Gi>Gi ^^^^ 

where ^'^ is the nominal Gibbs free energy of activation. Inserting Eq. 39 into Eq. 38, 
the forward and reverse kinetic rates are given by 

ki = vi exp(-/37*) 
A;_i = z/_i exp[-/3(7t - AG)] ^ ' 
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with forward and backward barrier crossing attempt frequencies ui and U-i. AG, 
AGj^^, ki, and k-i are shown in supplemental Fig. 1. From detailed balance, 
k-i/ki = exp(/3AG), and we find that the forward and backward barrier crossing 
attempt frequencies are identical, ui = U-i = u. 

The observed relaxation rate in any perturbation of a system governed by 
Eq. 36 is the eigenvalue = fci + k-i. Substituting the expressions for ki and 
k-i and rearranging, we obtain, = u e-x-p{—/3'y^)[l + exp(/3AG)]. The change 
in Gibbs free energy AG depends on the change in enthalpy AH and the change 
in entropy AS' through, AG = AH — TAS. Similarly, the change in the nominal 
activation free energy 7-^ depends on the nominal change in the enthalpy of activation 
5^ and the nominal change in the entropy of activation through, 'j^ = — Ta^. 
To simplify the derivation we ignore any temperature dependence in Ai7, AS, 6^^ 
and over the limited experimental range of temperatures. Making the substitution, 
P^t = _ pTai = + ai/R, we obtain 

r~ = u exp ( a ^r) exp(-/35*)[l + exp(/3AG)], (41) 



where z/adj = ^exp (cr-^/i?) is the forward entropy- adjusted (effective) barrier crossing 
attempt frequency. Taking the logarithm and defining = exp(/3AG) = k-i/ki = 
we obtain the equation for the curve of a reversible reaction in an Arrhenius plot 

ln(r-^) = In(^) + a^R - (36^ + ln(l + x). 

Recalling, 91n(l + x) = {1 + x)~^d{l + x)^ and observing, dx/df3 = xAH^ we obtain 
the desired expression for the slope in an Arrhenius plot from the relaxation rate of a 
reversible reaction 

In the forward irreversable limit {k-i <C fci), we obtain the classical result, dlB.{X)/d/3 = 
—S^. In the backward irreversable limit {k-i ^ /ci), dln{X)/d(3 — —S^ + AH. In a 
fully reversable reaction, {k-i = /ci), dln{X)/dp = —6^ + AH/2. In an Arrhenius anal- 
ysis that ignores reversibility, activation enthalpy is underestimated for endothermic 
reactions {AH > 0), and activation enthalpy is overestimated for exothermic reactions 
{AH < 0). 



Appendix B 

The equilibrium association constants Kq, i^i, and K2 for Ca^"^ binding to the system 
states (\00) , (\10), and (\11) are defined based on two assumptions, (i) The free 
energies of all ligand-unbound species Gq\s^s^ depend linearly on the chemical poten- 
tial of [Ca^^], /i = i?T ln[Ca^^], and the free energy of the protein complex Gq^^^^^: 
Go\siS2 — ^o\siS2 + (^ = !)• (ii) The free energies of the ligand bound species 
Gi\siS2 fixed: Gi^s-^^^ — G^^^^^^ (n = 0). For Ca^"^ binding to some configuration 
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of TnC-Tnl {\s1S2), Eq. 9 gives 





— 



0\S1S2 



The equilibrium association constant Kq for 



Ko = W+/W- = 




and 



Thus, i^o[Ca^^] is the pseudo-equihbrium constant that governs Ca^^ binding to fuhy 
deactivated TnC-Tnl (i.e. the (0\00) (1\00) transition). i^^i[Ca2+] and K2[Cd!^^] 
are the pseudo-equihbrium constants that govern Ca^^ binding to partiahy activated 
(0\10) and fully activated (0\11) TnC-Tnl. Ki and K2 are related to Kq by macroscopic 
detailed balance (Eq. 10): K0W1W-3 = K1W3W-1 and K1W2W-2 = K2W2W-2, 
which simplifies to Ki = K2. 
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Table 1: Microscopic parametrization of the Ca -signaling cycle. Parameter 
values are from a global fit of activation and deactivation stopped flow FRET 
data (c.f. Figs. 5 and 6) to the model of the signaling cycle (Fig 2, Eqs. 14 
and 15) (15 C). Reported values are mean ± precision. Precision is expressed as 
{/e, ue}^ where le and ue are the respective lower and upper 68% (1 cr) confidence 
estimates of the parameter, obtained from a 1-D adiabatic grid search. ^ fixed 
(see text). ^ poorly resolved, confidence estimates unavailable. ^ precision not 
determined because value depends on assumed ^ fixed to value measured 
from time-resolved FRET {dio ± a). 



Parameter 


global fit 


h is-') 


18.0 ± {14.4,21.9} 


k-i 


1.0" 




181 ±{145, 222} 




117 ±{93, 148} 


h 


1618'' 


k-3 


31,000« 


doo (A) 


21.4^ 


dio 


23.9 ± 2.6'^ 


dn 


32.2 ±{31.3, 34.5} 
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Table 2: Thermo- kinetic parametrization of the opening/closing transition of 
cardiac TnC-Tnl. For opening: Ai7, net enthalpy change; AC^, heat capacity 
change; AS^ net entropy change; S^^ enthalpy of activation; z^adj = z^exp {a^ /li)^ 
entropy- adjusted barrier crossing attempt frequency (z/, barrier crossing attempt 
frequency; cr^, entropy of activation.) For closing: enthalpy of activation; and 
^adj,R entropy-adjusted barrier crossing attempt frequency. 



Parameter 


Value 


A// (kJ/mol) 


-33.4 




10.6 




44.0 


(kJ/mol/K) 


-0.110 




-7.6 




1.8 X 10^ 


'^adj.R 


1.0 X 10^° 
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Figure Legends 

Figure 1. 

FRET assay for the structural kinetics of the CRS. The engineered FRET donor Trpl2 
(D), the FRET acceptor Cys51-labeled AEDANS (A), and the N-terminal portion of 
TnC (ribbon) are shown. The efficiency of photon transfer between the emission dipole 
of Trpl2 and the absorption dipole of AEDANS (shown) depends on the inter-probe 
distance (R), which changes upon hinge- like movement of helix B of TnC (q^b) during 
activation, (a) Ca^"^-bound TnC in the absence of Tnl (modified from PDB ID code 
1AP4). (b) The Ca^+-bound TnC with bound Tnl-R (brown helix) (modified from 
PDB ID code IJIE). 

Figure 2. 

The TnC-Tnl assembly represented as a Markov network (see text), (a) The parent 
network with system-states {so\siS2) as nodes and bi-directional transitions as edges, 
(b) Reduced models of activation and deactivation showing transition rates ki^i = 
1,2,3, Ca^^ affinity constants Kj^j = 0, 1,2, kinetic linkages between activation and 
deactiation, and degeneracy in k2 and k-2. Dominant (red) are minor (blue) pathways 
and common portions (green) are shown. Double arrows (<^) indicates transitions with 
the reversible binding of Ca^"^. For protein isomer izat ions, single arrows indicate the 
direction of net probability flux. 

Figure 3. 

Recovered FRET distance distribution. Background-corrected time-resolved donor- 
only and donor-acceptor decays were fit to a Gaussian distributed population of inter- 
probe distances with mean (bar) and standard deviation (shown as an error bar). 
Samples: starting Ca^^-free samples of TnC-Tnl complex or isolated TnC (bold) with 
addition of saturating Ca^^ or sequential addition of saturating Ca^^ then Tnl. 

Figure 4. 

Mechanistic models of TnC-Tnl activation consistent with the distance changes in 
Fig. 3. (a) Population-shift-stabilization model: TnC first undergoes an energetically 
unfavorable opening event (1), which is subsequently stabilized by Tnl-R binding (2). 
(b) Induced-fit model: Ca^^ binding produces an energetically favorable Ca^^-primed 
species that is not open (1). In the presence of Tnl, the Ca^^-primed complex under- 
goes concerted TnC opening and Tnl-R binding (2). The kinetic signatures of each 
model are shown, with relative arrow length indicating relative transition probabilites. 

Figure 5. 

Activation kinetics. FRET optical distance versus time, 15 C. (a) Ca^^-induced ac- 
tivation: change following rapid mixing of preformed binary TnC-Tnl complex with 
Ca^"^ solution (top dotted trace) or with buffer (mock injection, bottom trace). Empiri- 
cal double exponential fit (green): r^^ = 1554 s~-'^(48% amplitude), = 19 s~-'^(52%). 
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Global model-based fit (red) with derived parameters in Table 1. (b) Tnl-induced acti- 
vation: change following rapid mixing of Tnl into Ca^^ -pre-saturated TnC. Empirical 
single exponential fit (green): = 305 s~^. Global model-based fit (red) with derived 
parameters in Table 1. 

Figure 6. 

Deactivation kinetics after Ca^"^ removal, 15 C. FRET optical distance change after 
rapid mixing of Ca^^-saturated TnC-Tnl with buffer containing Ca^^ chelator, EGTA. 
Empirical single exponential fit (green, obscured by red trace): = 125 s~^. Global 
model-based fit (red) with derived parameters in Table 1. Inset: (Eq. 22), the 
normalized adiabatic projection of the global hyper-surface on a model parameter 
(labeled) axis. 

Figure 7. 

System response to a change in [Ca^^]. (a) Ca^^ titration, mean FRET distance vs. 
pCa (-log([Ca^^])). Empirical fit to the Hill equation (green) (n = 1.18; pCa^o = 5.82). 
Global fit to the model (Eqs. 2, 5, 6) (red) (/3A/i' = 22.9 kJ/mol). (b) The experiment- 
resolved BFEL (see text, T = 15 C). Ca^^-unbound system-states (green). Ca^^-bound 
system-states (blue). Population density (spheres, radius represents magnitude) at the 
extremes of activation. pCa^Q (labeled, ►). 

Figure 8. 

Arrhenius analysis of the induced fit opening step of activation. The rate of Tnl- 
induced activation (circles), r"^ = A:2 + at T = {3.7,10,15,16,20} C (c.f. Fig- 
ure 5b) were fit to Eq. 26 (solid line) to recover the enthalpy of activation and 
the entropy- adjusted barrier crossing attempt frequency z/adj foi" the opening transition 
(Table 2). Also shown (squares) are the k2{T) (filled) and k-2{T) (open) recovered 
from global analysis of relaxation data at T = {3.7, 10, 15, 16} C (c.f. Figures 5 and 6) 
and calculated k2 and k-2 (dashed lines) from the recovered z/adj? ^adj R' T''"? 7r- 

Figure 9. 

Thermodynamic analysis of the induced fit opening step of activation. The temperature- 
dependent equilibrium constant, Ka = ^^2/^-2, (circles) from k2{T) and k-2{T) (Fig. 8) 
were fit using Eq 25 with temperature-independent Ai7 and AS' (solid line) and 
temperature-dependent Ai7 and AS (dashed line). Recovered Ai7, AS, and Cp 
are given in Table 2. Inset: van't Hoff plot showing extrapolated results from fit 
to temperature- independent Ai7 and AS. The y-intercept and slope provide —AS/R 
and Ai7, respectively. 

Scheme 1. 

Coarse-grained model of an allosteric protein. 



dynamic equilibrium 37 

Figures and Schemes 



dynamic equilibrium 




Figure 1 



dynamic equilibrium 



1\01 



1\11 



1\00- 



1\10 



1\00- 



1\11 



0\01 



-0\11 



^l\10k. 



0\00- 



0\10 



0\00- 



■0\10 0\00- 



s 

i\io ■ 

K 



1\11 



'0\11 

V 



■o\io K 



activation 



deactivation 



Figure 2: 



dynamic equilibrium 



^ 32 



c 

'-t—f 

-I— • 
CO 

b 


o 
c 

w 

b 
I- 

LU 
DC 



30 
28 
26 ^ 
24 
22 



20 



22.3 
(1.6) 



28.9 
(2.9) 



21.7 
(2.6) 



23.9 
(2.6) 



28.3 
(2.8) 



Figure 3: 



dynamic equilibrium 




dynamic equilibrium 




Figure 5: 



dynamic equilibrium 




10 20 30 40 50 

Time (msec) 



Figure 6: 



dynamic equilibrium 



(a) 

30 1 




(\11) 8 

Figure 7: 



dynamic equilibrium 




Figure 8: 



dynamic equilibrium 




0.1 0.2 0.3 0.4 

1/RT (kJ/mol)-^ 



0.415 0.420 0.425 0.430 
1/RT (kJ/mol)-^ 



0.435 



Figure 9: 



dynamic equilibrium 



47 



■□ ^ HO 

Cn fco 

mm 

Scheme 1: 



dynamic equilibrium 48 



Supplemental Online Data 
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Supplemental Figure 1: Energy surface of a reversible noise-induced transition 
in a bi-metastable potential. 




Supplemental Figure 2: Time resolved decays of Trpl2 in 
TnC(12W/C51±AEDANS) showing the effect of Tnl addition (48.8 ps/channel). 
TnC in 160 fiM Ca^^: donor-only (black), donor-accepter (green). TnC in 
160 fiM Ca^^ + 2-fold excess of Tnl: donor-only (red). Donor-acceptor (blue). 
Impulse response (black line). 



